Unit- II 


Linear Systems- Let us consider a system of first order differential equations of the form 

( 1 ) 


£ =FU.y ) 

dt 


< Y = G(x,y) 
dt 

Where t is an independent variable. And x & y are dependent variables. 

The system (1) is called a linear system if both F(x, y) and G(x, y) in x and y. 


Also system (1) can be written as 


dx 

— = a l (t)x + b l (t)y + f l (t) 
dt 


^j~ = a 2 ( t)x + b 2 ( t)y + f 2 (t) 
dt 


( 2 ) 


Where a ( .(f), b^t) and a { {t) Vi = 1,2 are continuous functions on \a,b\. 

Homogeneous and Non-Homogeneous Linear Systems- The system (2) is called a 
homogeneous linear system, if both /' (7)and f 2 (t) are identically zero and if both f\ (7) and 
f 2 (t) are not equal to zero, then the system (2) is called a non-homogeneous linear system. 

x = x(t) r , 

Solution- A pair of functions defined on [a,/?Jis said to be a solution of (2) if it satisfies 

y = y(0 

( 2 ). 


Example- 


dx 

— = 4x - y 
dt 

— = 2x+y 
dt 


..A 

.B 


(3) 


. . dx . . , d X — dx . _ . ^nd i i • rr • i 

From A, y = 4x -putting in B we obtain —— - 5-f 6x = 0 is a 2 order differential 

dt dt dt 

x — o 

equation. The auxiliary equation is m - 5m + 6 = 0 => m - 2, 3 so putting x = e in A, we 


x = e 


obtain y = 2e ir again putting x = e 3 ' in A, we obtain y = e i '. Therefore the solutions of (3) are 


2 1 3 1 

x=e x=e 

and 

y-2e y - e 


( 4 ) 





Theorem ! If t Q is any point of [a,Z?]and x 0 & y 0 are any two numbers, then the system (2) has a 

x = x(t ) x(t 0 ) = x Q 

unique solution with 

y = y(t) y(t 0 ) = y 0 


dx 

Theorem-2 If the homogeneous system 

dy 


dt 


a x (t)x + b l ( t)y 
a 2 (t)x + b 2 (t)y 


x = x 1 (t) X = x 2 (t) 
has two solutions and 

y = yi(t) y = y 2 ( 0 


(5) 


( 6 ) 


r n X = C.X,(t) +C^xJt) 

on \a,b\. Then 11 22 

y = c 1 y l (t) + c 2 y 2 (t) 


is also a solution of (5) on [a,Z?]for any two constants c { and 


(7) 


x = x l (t) x = x 1 (t) 

Theorem-3 If the two solutions and 

y = yi(t) y = y 2 ( 0 


(6) of the homogeneous system (5) 


r x = CjXj (t) + c 2 x 2 (t) 

have a wronskian W(l) that does not vanish on [a,b \, then (7) is a general 

y = c 1 y l ( t) + c 2 y 2 (t) 

solution of homogeneous system (5) on [a,b\. 


Note- The wronskian W(t) of the solutions (4) is 

3 1 2 1 


W(t) = 


e 

5 1 


e 

2e 2 ' 


= e 

Theorem -4 The wronskian W(t) of two solutions (6) of homogeneous system (5) is either 
identically zero or nowhere zero on [a,Z?]i.e 

W(t) - 0(linearly dependent) or W(!) ^ 0 (linearly independent). 


The wronskian W(t) satisfies the differential equation, 
between the limits 0 to t we obtain 


dW 

dt 


[a 1 (t) + Z? 2 (t)]Vk and on integrating 


W ( t ) = cej [aj (, t) + b 2 ( t)\it 
o 





X = X, (t ) X = x 2 (t) 

Theorem -5 If the two solutions and of homogeneous system (5) are linearly 

y = yi(t) y = y 2 (0 

x = .r p (0 

independent on [a, b\ and if is any particular solution of non-homogeneous system (2) 

V = .'•„(/) 

jc = c 1 ac 1 (0 + c 2 ac 2 (0 + ^ / ,(0 

on [a, b\ , then is a general solution of of non-homogeneous system 

y = c l y 1 (t) + c 2 y 2 (t) + y p (t ) 

(2) on [a, Z?]. 


Example- Show that 


At - 2 1 

x = e x = e 

and are the solutions of the homogeneous system 

y = e 4t y = —e~ 2t 


dx_ 

dt ^ x = x (0 

and find the particular solution of the given system for which jc( 0) = 5 and 

dy , , y = y(t) 

— = 3x+ y 

dt 

>’( 0 ) = 1 . 


Solution- Let 


— = x + 3 v 

dt 

— = 'ix+y 

dt 


At - 2 1 

x = e x = e 

First, we show that each of the pair and ^ satisfy the system (1). In order to 

y = e 4t y = -e~ 2t 

X = CjXj (t ) + c 2 x 2 (t ) 

determines a particular solution of (1), let us consider (2) be a particular 

y = c 1 y l (t) + c 2 y 2 (t) 

solution of (1), where the constants c x and c 2 are to be determined. Putting the values of 
x 1 (t) = e 4> , x 2 (t) = e~ 21 , y 1 (t) = e 4 ' and y 2 (t) = -e~ 2t in (2) and using the given conditions 

x(0) = 5 and y(0) = 1, we obtain c 1 = 3 and c 2 = 2. 


x = 3e 4t + 2e~ 2t 

Therefore is a particular solution. 

y = 3e —2e 


x = 3t - 2 

Example Show that is a particular solution of the non-homogeneous system 

y = -21 + 3 



dx 

— = x+2y + t-l 
dt 

< 

— = 3x+ 2y-5t- 2 
dt 


and write the general solution of this system. 


Hint- Let 


dx 

— = x + 2y + t-l 
dt 

— = 3x + 2 y - 5f - 2 
dt 


x = 3t-2 

Now will be a particular solution of the non-homogeneous system (1) if it satisfies 

y = -21 + 3 

the system (1). In order to find a general solution of system (1), we have to find a solution 

dx 

— = x + 2 y 

corresponding homogeneous system (2) to system (1) as similar in example in 

— = 3x+2y 
dt 

equation (3). 


Answer- 


x = 2c x e 4 ' + c 2 e ’ + 3t - 2 
y = 3c t e 4 ' -c 2 e~' — 2t-\-3 


Homogeneous Linear Systems with Constant Coefficients- Let us consider a homogeneous 


linear system with constant coefficients 


dx 

— = a { x + b { y 
dt 

dy , 

— = a 2 x + b 2 y 

dt 


x = Ae mt 

Where a,,b,, a-, and b n are constants. Suppose (2) 

i i z z a a mt 

y- Be 

(where A, B and m are to be determined) be a solution of the system (1), then it satisfies (1) so 
Ame mt =(a l A + b i B)e" ,t 
Bme mt = {a 2 A + b 2 B)e mt 


(a l - in)A + b x B = 0 
a-,A + (b 2 —m)B = 0 


( 3 ) 



is a system of equations of the form ax — 0 has a trivial solution x = 0, if A-B- 0 so for a 
nontrivial solution x A 0 of (3), we have a - Oi.e 


b 2 - m 


= 0 , on expanding we obtain a quadratic equation in m 


m — (<a, + b 2 )m + (a l b 2 - a 2 b x ) - 0 

gives two values of m say m l and m 2 . Now the following three cases arise 


Case-1 If m 1 and m 2 are real and distinct, then corresponding to m l , we find the values of A and 

x = A.e m ' t 

T~» _ A _1 n 1_ _ A.' _ _ A. 1- _ _ J. __ • 1 O!_! 1_ 


B say A, and by equation (3), so the first nontrivial solution is 


y = B { e’ 


Similarly 


corresponding to m 2 , we find the another nontrivial solution 


A m? t 
2 e - 

r> T 

y = B 2 e 


Therefore the general solution is 


x = c l (A l e m ' r ) + c 2 (A 2 e m2 ') 
y = c 1 {B l e m ' t ) + c 2 (B 2 e m > t ) 


Example- Find the general solution of the system of equations 


■ = x + y 


— = 4x~2y 
dy 


Solution- Let 


— = x + y 
dt 

— = 4x~2y 
dy 


On comparing a { = Y,b x = l,a 2 = 4and b 2 =— 2, the auxiliary equation is m 2 +m -6 = Ogives 
m = -3, 2 


Where A and B satisfy 


(1 - in) A + B = 0 
4A + (-2 - m)B = 0 


When m - -3, then by (2) we get A = 1, B = -4 and the first nontrivial solution is 


y = -4A 



Similarly for m = 2 , then by (2) we get A = 1, B = 1 and the another nontrivial solution is 


y = 4e 1 


Therefore the general solution is 


x = c x e 31 +c 2 e^‘ 
y = -4c t e 3 ' +c 2 e 2t 


Example- Find the general solution of the system 


dx . , 

— = -3 + 4v 
dt 

^=-2j + 3y 
dy 


Answer 


x = 2 c { e ' +c 2 e' 
y = c x e~' +c 2 e‘ 


Case-2 If m l and m 2 are conjugate complex numbers of the form a + ib, where a and b are real 

x= Ay a+ib)t 

numbers with b ^ 0, then we consider two linearly independent solutions (1) and 

y = By a+ib)t 

x = A*e'' a ~ ib)t 

( ^ , where A~ = A t + iA 2 , B x = B x + iB 1 ,A* 2 =A X - iA 2 and B* = B x - iB 2 resp. Putting 


y = 

the values of A, and B x in (1), we have 

x = (Aj + iA 2 )e at (cosbt + isinbt) 
y = (5, + iB 2 )e al (cos bt + i sin bt) 


x = e"'[(Aj cos bt- A 2 sinZtf) + /(A 1 sinZ?t + A 2 cosZ?r)] 
y = e a, [(B l cos bt- B 2 sin bt) + i(B x sin bt + B 2 cosZ?t)] 


Equating real and imaginary parts, we obtain two linearly independent solutions say 


x = e‘" (A, cos bt-A 2 sin bt) 
y = e at (B { cos bt-B 2 sin bt) 


(3) and 


x = e"'(A, sin bt-A 2 cos bt) 
y = e at (B { sin bt-B 2 cos bt) 


Therefore the general solution is 

x = e a, \c x (A x cos bt-A 2 sin bt) + c 2 (A x sin bt + A 2 cosZtf)] 
y = e a, \c x {B x cos bt — B 2 sin bt) + c 2 (B l sin bt + B 2 cosZ?t)] 



Example- 


dx 

— = 4x-2y 
dt 

— = 5x + 2y 
dt 


dx 

- = 4x-2y 

Hint- dt 

— = 5x + 2y 
dt 


The auxiliary equation is m 2 — 6m + 18 = Ogives m = 3±3i, taking a nontrivial solution 
x = (A, +/A 2 )e 3 '(cos3t + /sin30 

(2) of (1) , where A l ,B l , A 2 and B ,are to be determined. For 

y = (B l +iB 2 )e 3 ‘(cos3t + isin3t) 

this (2) satisfies (1) and equating the coefficients of cos 3t and sin3ion both sides. 


Answer- 


x = e 3, (2c 1 cos3t + 2c 2 sin3t) 
y = e 3 ' [cj (cos 3t + 3 sin 3 1) + c 2 (sin 3t — 3 cos 3t)] 


Case -3 If m { = m 2 = m are equal roots then we should have only one linearly solution 

x = Ae mt nd x = Ate"" 

and the 2 linearly independent solution will be of the form . But actually, 

d mt r>. mt 

y- Be y — Bte 

we consider the 2 nd linearly independent solution 


x = (A, + A 2 t)e" 
y = (B, + Bj)e' 


t , where A,B,A l ,A 1 ,B l and B 2 are to be determined. 


x = c l Ae mt + c 2 (Aj + A 2 t)e" 

Therefore the general solution is 

y = CjBc"" + c 2 ( B , + B 2 t)e" 


Example- Find the general solution of the system 


— = 3jr-4y 

dt 



Solution- Let 


( 1 ) 


dx 

dt 


3x-4 y 



y 


The auxiliary equation is 
m 2 —2m + l = 0 
m = 1,1 


x = Ae 
Let 

y = Be 

be a solution of (1), where A and B satisfy 


( 2 ) 


2A-4B = 0 


A-2B = 0 


gives A = 2, B — 1 ,so 


x = 2e 
y = e 


(3) 


be a first linearly independent solution of (1). We consider the second linearly independent 

x = (A| + A-,t)g 


solution of (1) of the form 


so it satisfies (1) 


y = (B 1 + B 2 t)e‘ 


(4) 


(2A, - A 2 - 45,) + (2A, - 4B 2 )t = 0 + 0t 

on equating both sides we have 

(Aj -2B i -B 2 ) + (A 2 -2B 2 )t = 0 + 0i 


2A.-A,-4fi,=0 A-2B.-B, =0 

1 2 1 and 1 1 2 (5) 

24-4^2=0 A 2 - 2B 2 = 0 

On solving the equations in (5), we obtain A l = 1, B { = 0, A 2 = 2 & S 2 = 1 

The another linearly independent solution is 

x = (l + 2t)e‘ 
y = te' 


Therefore the general solution is 



x = 2c,e r + c 2 (1 + 2 t)e' 
y = c { e' +c 2 te r 

Example- Find the general solution of the system 


dx 

dt 


= 5x + 4 y 


dy 

dt 


= —x + y 


x = -2c { e 3 ' +c 2 (l + 2t)e 3 ' 

Answer 

y = c x e 3t -c 2 te3t 

Non-Linear Systems: Volterra’s Prey- Predator Equations- 

Everyone knows that there is a constant struggle for survival among different species of animals 
living in the same environment. One kind of animal survives by eating another and a second by 

For an example of this universal conflict between the predator and its prey, let us imagine an 
island inhabited by foxes and rabbits. The foxes eat rabbits and the rabbits eat clovers. Let us 
assume that there is so much clovers then the rabbits have an ample supply of food. When the 
rabbits are abundant, then the foxes flourish and their population grows. When the foxes become 
too numerous and eat too many rabbits, then they enter into a period of famine and their 
population begins to decline. As the foxes decrease, then the rabbits become relatively safe and 
their population starts to increase again. Thus we have an endless repeated cycle of the increase 
and decrease in two species of animals and the fluctuations in two species are given by the 
following figure 



If x and y are the number of rabbits and foxes at any time t, then in the presence of an unlimited 
supply of clovers, 




dx 

The rate of change of rabbits is — = ax, a > 0, after some encounter between the rabbits and 

dt 


dx 

foxes the rate of change of rabbits is — = ax-bxy, a,b> 0 

dt 


( 1 ) 


In the absence of rabbits the foxes die and the rate of change of foxes is — = -cy, c > 0 and 

dt 

after some encounter of foxes with rabbits their population grows and the rate of change of foxes 
become 


— = -cy + dxy, c,d > 0 
dt 


( 2 ) 


These two equations are called the volterra’s prey-predator equations. 
For the solution of these equations, we divide (2) by (1) 
dy 

dt _ ~y{c~dx ) 
dx x(a-by) 
dt 


Or 

dy y{c - dx) 
dx x(a-by) 

on separating the variables, we have 

(c - dx)dx + (a-by)dy _ 
v y 


f —d dx+ [- b 

•M v ) J\y 


kv=o 


On integrating, we have 

c log x + a log y = dx + by + log K 

or x c y a = Ke (dx+hy) 


(3) 


( 4 ) 


In order to determine K putting x(t 0 ) = x 0 , y(t 0 ) = y 0 in (4) so 



Therefore the solution of volterra’s prey- predator equations is 


x y 


c ~(dx 0 +by Q 


v 0 


y 0 e 


j^(dx+by) 


Non-Linear Equations- Let us consider the motion of a pendulum consisting a bob of mass m 
attached to one end of a light rod of length a. If the bob is pulled to one side through an angle a 
and and then released, let 6 be the position of the bob after time t s.t. AQ= s, thgn by the 


o 



principle of conservation of energy, 

Gain in kinetic energy =Loss in potential energy 


1 2 

—mv 

2 



= mg (a cos 6- a cos a) 


ga(cosQ-cosoc) 


( 1 ) 


. - „ ds dO 

Also s -aa, so v = — = a —, 
dt dt 


putting in (1) 


1 2 (cW\ 2 
— a" — 

2 \dt) 


ag{ cos#-cos a) 


1 (d0\ 

— a - 

2 k dt 


2 


g (cos d?-cos a) 


Differentiating w. r. to t 















f 1 ^ 

f r. d6\ 

f d 2 e) 

r 

—a 

2— 


= g 

U J 

v dt j 

v dt 1 j 

V 


• Q dO 

-sin 0 -0 

dt 


d 2 0 g . a 
— r = --sin0 
dt~ a 

Replacing 6 by x, so 


d 2 x g ■ n 
dt a 

is a non-linear differential equation of first order . 


( 1 ) 


d x p 

If x is small, thensinx = x, so it becomes linear—— + — x = 0, if the damping (or resistance) 

dt a 

force is proportional to velocity, then the equation of motion is 


d 1 x ( c x 


dt 2 


■ + 




dx g . _ 

-h— sinx = 0 

dt a 


is a non-linear differential equation of 2 nd order. 


( 2 ) 


ai d 2 x . 2 , dx 

Also— —+u(x -1)-l-x = 0 

dt 2 dt 


(3) 


is a vander pol equation. 


Now, we consider a 2 nd order non-linear differential equations of the form 


d x 

h ^ 2 


= f 


dx 

x, — 

dt 


(4) 


Autonomous System and Phase Plane- Suppose a particle of unit mass moves on the x axis 


and / 


dx 

x ’ T 
V dt j 


( dx 

is the force acting on it, then the values of pair x{position), — ( velocity) 

v dt j 

dx 


are 


called the phase of the system at each instant and the plane containing x and — is called the 
phase plane. 

Introducing the new variable y = —, then the 2 nd order non-linear differential system (4) is 

dt 

equivalent to the system 



dx 


y 


dt 

d 2 x d 


i.e 


dt 2 dt 


= -(y) = f(x, y) 


dx 

yr y 

^-=f(x, y ) 
dt 


(5) 


A system of the form — = f{x(t), y(t)) is called a non-autonomous differential system and a 

dt 

system — = f(x, y) in which the RHS does not contain the independent variable is called a 
dt 

autonomous differential system. 


The functions x(t) and y(t) (where t is a parameter) are the solutions of (5) and define a curve in 


the x-y plane, which is also called the phase plane because 


^ dx ^ 

x,y = — 

\ at j 


Now in general, we consider a system of the form 


d -d=Fix,y) 

dt 

~ = G(x,y) 
dt 


( 6 ) 


Where F(x,y) and G(x,y) are continuous functions of x & y and have continuous first partial 
derivatives in the phase plane. 

Note- If t 0 is any number and (x 0 , y 0 )is any point in the phase plane the there exists a unique 

x = x(t) x(t n ) = x n 

solution of (6) with 

y = y(t) y(t 0 ) = y 0 

X = x{t) 

Path- If both the functions x(t) and y(t) are not constant functions, then the solution of 

y = y(t) 

(6) is a curve in the phase plane and it is also called a path of the system. 


x = x(t) x = x(t) 

Note- If is a solution of (6) then is also a solution of (6) for any constant c. 

y = y(t) y = y(t ) 

Thus each path is represented by many solutions which differ from one another only by a 
translation of parameter. 


Critical Points- The points (jc 0 , y 0 ) at which both the functions F and G vanish 



i.e 


F(x 0 , y 0 ) = Oand G(x 0 , y 0 ) = 0are called the critical points of the system 


— = F(x, y ) 

at , 

and 

C y = G(x,y ) 
dt 


JV - ^VQ 

at such a point a unique constant solution exists and it is not a path. Hence no path goes 

y = y 0 

through a critical point. 

Isolated Critical Point- A critical point (^c 0 , _y 0 ) is said to be an isolated critical point if there 
exists a circle centered on (jc 0 , y 0 ) that contains no other critical point. 

Note- The followings are needed to describe a phase portrait of two dimensional fluid motion, 

1. The critical points. 

2. The arrangement of paths near critical points. 

3. The stability of critical points. 

4. Closed paths. 

Example- Describe the phase portrait of the following 


dx 

dx 

dx 

— = 0 

— = x 

— 

dt 

dt 

dt 

|£ 

ll 

o 

2 -ffe = o 

dy 

dt 

dt 

dt 


Solution- 1.For critical points putting each of F and G equal to zero so F(x, y)=0 and G(x, y)=0 
for all values of x and y. Therefore each point of a phase plane is a critical point. For paths either 

dy 

-G( X v) -dL- G(x,y) 

and ~^x,y) t 


dx 

— — F(jc y) 
we integrate separately ^ ’ - ’ 


dt 


or we integrate dx 
dt 


F(x,y) 


^ = 0& ^ = 0 

So on integrating j t j t we have 


x = Cj & y = c 2 


. Hence no path exists(because both 


x and y are constant functions). 


Answer 2- every point on y axis is a critical point and paths are horizontal half linesdirected out 
to the left and right from y axis. 


Answer 3- point (0, 0) is the only critical point, and paths are half lines of all possible slopes 
directed in toward the origen. 


Example- Find the critical points of the following 



1 - 


2 - 


d 2 x dx , 2 „ , 

—- +-(x 3 + x 2 -2x) = 0 

dt dt 

dx 2 z r 
— = y -5x + 6 

dt 
dy 


dt 


= x-y 


dx . dx dx 


dt " dt 2 dt 


Solution- 1 putting = y in - (x 3 + x 2 - 2x) = 0, so it becomes 

dx 


dt 


= y 


— = (x 3 + x 2 - 2x) - y 
dt 


( 1 ) 


For critical points putting y = 0 & (x 3 + x 2 - 2x) - y = 0, we have y = 0 & x(x 2 + x - 2) = 0 
x = 0,1, - 2 

So the critical points are (0,0), (1, 0) & (-2, 0). 


Answer-2 (2, 2) & (3, 3). 




Unit- III 


Critical Points & Stability for Linear Systems- Let us consider an autonomous linear system 
with constant coefficients 


dx 

— = F(x,y ) 
dt 

-y = G(x,y ) 
dt 


( 1 ) 


Then (1) can be written as 


dx . 

— = a x x + b { y 
dt 

dy 

— = a 2 x + b 2 y 
dt 


( 2 ) 


x = Ae mt 

has a critical point (0,0). Suppose be the solution of (2) and if 

y = Be mt 

auxiliary equation of system (2) is 


h x 

b. 


^ 0, then the 


in - (a, + b 2 )m + ( a x b 2 - a 2 b { ) = 0 (3) 

and it is a quadratic equation in m. The nature of a critical point of (2) is determined by the 
nature of roots m 1 & m 2 of (3) which can be classified in two categories such as major cases and 
minor cases. 


Major Cases: Case 1- The roots m 1 & m 2 are real, distinct and of same sign, then the critical 
point is a node i.e 

m 2 ^ m 2 

- - {Node & asymptotically stable) 

+ + {Node & unstable) 


Case 2- The roots m x & m 2 are real, distinct and of opposite signs, then the critical point is a 
saddle point i.e 

m x -£■ m 2 

- + {saddle ^unstable) 


Case 3- The roots m l & m, are conjugate complex numbers, but not purely imaginary, then the 
critical point is spiral i.e 


{ 16 } 



- a + ib(spiral & asymptotically stable) 
a - ib (spiral & unstable) 


m x & m 2 = 


and — = ± 
dt 


Borderline cases: Case 4- The roots m { & m 2 are real and equal, then the critical point is a node 
i.e m { = m 2 = m < 0 (node & asymptotically stable) and if 

m ] = m-, = m > 0 (node & un stable). 


Case 5- The roots m x & m 2 are purely imaginary, then the critical point is a centre and stable, but 
not asymptotically stable. 

Example- For each of the given linear system , find (1) Critical points (2) find the general 
solution (3) find the differential equation of paths (4) solve the equation of paths (5) sketch a 
few of the paths (6) discuss the stability of critical points. 


a- 


dx 

dt 

dy 

dt 


= -x 


b- 


dx 

— = 4 y 

dt 
dy 


— = —2 y — = -x 

dt 


x — c { e dy 2y 

Solution: a- Flere (1) the critical point is (0, 0), (2) solution is , (3) — = —, 

y = c 2 e~ 2 ' dx x 


(4)y = cx 2 , (5) paths are exponential curves, (6) since both the limits 
lim M ^ x(t) and lim,^ y(t) exist, then the critical point is asymptotically stable. 


Note- If lim^^both x(/)and \'(/)exist, then the critical point is asymptotically stable and if 
lim Moo one of these x(i)and y(t) exist, then the critical point is unstable. 


x = 2c, cos 2 1 + 2c 2 sin 2 1 


Solution b- (1) (0,0), (2) 

y = -c l sin 2 1 + c 2 cos 2 1 
stable but not asymptotically stable. 


(3) £= * (4) 

dx 4 y 4c 


2 2 

+—— = l. 

c 2 


(5) 


Note- Stable and unstable critical points: A critical point is said to be stable, if for each 
positive number R there exists a positive number r < R such that every path which is inside the 
circle x + v = r for some t — t Q remains inside the circle x +y = R~ for some t>t 0 . A 
critical point is said to be asymptotically stable, if it is stable and there exists a circle 



x 2 + y 2 = r 0 2 such that every path which is inside this circle for some t — t 0 , approaches the origin 
as t —» °° and if the critical point is not stable, then it is called unstable. 



dx , 

— = a { x + b Y y 

Theorem- The critical point (0,0) of the linear system is stable if and only if the 

dy , 

— = a 2 x + b 2 y 
dt 

roots of the auxiliary of the given system have non positive real parts, and it is asymptotically 
stable if and only if both roots have negative real parts. 

Proof- The auxiliary equation of the given system can be written as 


(m - m l )(m - m 2 ) = m 2 + pm + q — 0 

Where p = ) & q = m l m 2 , 

— pi- 

then by (1), we have m l , m 2 = - 


( 1 ) 

let us consider p and q axes and excluding the case q — 0, 

4p 2 ~ 4< ? 

2 


Now the three cases arise 


4p 2 -4 q =0, + 


or - 


If p 1 - 4q — 0, i.e on the parabola p 1 - 4q — 0 ,them the roots & m 1 are real and equal and the 
critical points are node. 

If p 2 -Aq <0, then the roots m, & in 2 are complex conjugate numbers, the critical points are 
spirals and centers. 









q 


Stable 



If p 1 - 4q > 0, then the roots m ] & m 2 are real, distinct and the critical points are saddle points. If 
p 2 - 4q > 0 and q< 0, then the roots m l & m 2 are real and distinct, then the critical points are 
spirals. 

dx , 

— = a x x + b { y 

Theorem- The critical point (0, 0) of the linear system is asymptotically stable if 

dx . 

— = a^x + b 2 y 
dt 

and only if the coefficients p = -( a { + b 1 ) & q = ( a l b 2 - a 1 b l ) of the auxiliary equation are both 
positive. 


Example- Determine the nature and stability properties of the critical point (0,0) for each of the 
following linear autonomous systems 


dx 

— = —x — 2y 
dt 

d ?=4x-5y 

dt 


2 . 


dx „ . 

— = -3x + 4y 
dt 

d >=-2x + 3y 
dt 


3. 


*=4x-3y 

dt 

& = Sx-6y 
dt 


{ » } 


1 . 








dx 

Solution 1-Let 

dy 


dt 


—x — 2 y 
4x-5 y 


— x — 2y = 0 

For critical points solving , we have (x, y) = (0,0) 

4x - 5 v = 0 


The auxiliary equation is m 2 - (-1 - 5 )m + (5 + 8) = 0 i.e m 2 + 6m + 13 = 0 


i.e m = —3 ± 2i , since m ] & m 2 are complex conjugate numbers a±ib & a < 0, then the critical 
point are spirals and asymptotically stable. 

Answer 2- Unstable saddle point 

Answer 3- The critical point is not isolated. 

Stability by Liapunov’s Direct Method- Liapunov’s direct method is used for studying the 

dx 

— = F(x,y) 

stability problems of a linear autonomous system (1) with an isolated critical 

^- = G(x,y) 
dt 

point (0,0). Let C = [x(t), y(t)]be any path of (1) and E(x,y) be any function which is 
continuous and has continuous first partial derivatives in a region containing this path. If a point 
moves along this path with x = x(t) & y = y(t) , then £(x, y) can be regarded as a function of t 
and the rate of change of £(x, y) with respect to t is 


dE _ dE dx ^ dE dy 
dt dx dt dy dt 

( 2 ) 

dx dy 

is called a Liapunov’s direct method. 

Positive Definite- A function E(x,y) is said to be positive definite, if 

£(0,0) = 0 & E(x, y)> 0 V (x, y) + (0,0). 

Negative Definite- A function E(x, y)is said to be negative definite, if 

E( 0,0) = 0 & E(x, y)< 0 V (x, y) * (0,0). 


Positive Semi Definite- A function E(x,y) is said to be positive semi definite, if 
£(0,0) = 0 & £(x, y) > 0 V (x, y) * (0,0). 



Negative Semi Definite- A function E(x, y)is said to be negative semi definite, if 
E( 0,0) = 0 & E(x, y)< 0 V (x, y) * (0,0). 


Note- Let us consider a positive definite function E(x, y ) of the form 

E(x, y ) = ax 2 '" +by 2n ,where a &b are positive constants and m & n are positive integers. 


Liapunov’s Function- A positive definite function E(x,y) is said to be a Liapunov function if 
dE dE dE . . . , ... . . dE . . 

— = — E H- (j is negative semi definite i.e — < 0, so it is a decreasing function. 

dt dx by dt 


Theorem- The function E(x, y) - ax 2 +bxy + cy 2 is positive definite, if and only if 
a >0 & b 2 - 4 ac < 0, and is negative definite if and only if a <0 & b 2 - 4 ac < 0 and is neither if 
a>0 & b 2 - 4 ac > 0 or a <0 & b 2 - 4 ac > 0 


Example- Determine whether each of the following function is positive definite, negative 
definite or neither: 

A- x 2 -xy-y 2 B- 2x 2 -3xy + 3y 2 C — 2x 2 +3xy-y 2 D— x 2 -4xy-5y 2 

Solution A- Here E(x, y) — x 2 - xy - y 2 comparing a = 1> 0 &b 2 -4ac = 5 > 0so it is neither. 

Answer B- Here E(x,y) = 2x 2 -3xy + 3y 2 comparing <2 = 2 > 0 &b 2 -4ac = -15 < 0, so it 
positive definite. 

Answer C- Here E(x,y) = -2x 2 +3xy - y 2 on comparing a = -2 < 0 &b 2 -4ac = 1 > 0, so it is 
neither. 


Answer D- Here E(x, y) = -x 2 - 4 xy - 5y 2 on comparing a = -l<0&Zr - 4 ac = -4 < 0, so it is 
negative definite. 


Example - Show that (0, 0) is an asymptotically stable critical point for each of the following 
systems: 


A- 


dx 3 

— = -ox - y 
dt 

dy 5 r, 3 

— = x - 2 v 
dt 


B- 


dx _ 3 

— = - 2 x + xy 
dt 

dy 2 2 3 

— = -x y -y 

dt 
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( 1 ) 


dx 

Solution- Let 

dy 


dt 


-3x 3 - y 
x 5 -2y 3 


On comparing F(x, y) - -3x 3 - y & G(x, y) = x 5 -2y 3 , Now we consider 
E{x,y) = ax 2 '" +by 2n where a,b,m&n are to be determined. The critical point (0,0) will be 


dE dE dE . . 

stable, it — = —F -1-G, is negative semi definite i.e 

dt dx dy 


— = 2amx lm ~ 1 (-3x 3 -y) + 2bny 2n ~\x 5 -2y 3 ) 
dt 


— = -6 amx 2m+2 - 2amx m 1 y + 2bnx 5 y 2n 1 -4bny ln+2 will be negative semi definite if 
dt 

- 2 amx 2m ~ l y + 2 bnx 5 y 2n_1 = 0 i.e 2amx 2m ~ l y = 2bnx 5 y 2 " _1 

equating the powers of x & y and also the coefficients of xy on both sides, we have 
2m — 1 = 5 => m = 3 & 2/7 -1 = 1 => n = 1 and 6a = 2b => a = 1 & b - 3 


So E(x,y) = x +3y . 

Similarly 2 nd part can be solved 
Answer 2- E(x, y) = x 2 + y 2 

Example- Show that (0, 0) is an unstable critical point for the system 


dx _ ■> 

— = 2xy + x 
dt 

dy 2 , 5 

— = -x + y 
dt 


Solution- Let us assume contradictory that (0, 0) is a stable critical point so 


— = 4 amx 2m y + 2amx 2m+2 - 2bnx 2 y 2n 1 + 2bny 2,1+4 will be negative semi definite if 
dt 


4amx 2m y = 2bnx 2 y ln 1 => m = n = \ & a = \ &b = 2, for these values of m,n,a&b , the value of 



— > 0 Hence the critical point origin is unstable. 
dt 


Green Function- Let us consider an nth order linear homogeneous differential equation 


L(y) = 0 (1) 

where L is a differential operator given by 

L (y) = P 0 (x)-^—+ p l (x)-^—+ +P n (x) (2) 

dx dx 

where p 0 (x), p { (x), p n (x) are continuous functions on \u, Z?]and the boundary 

conditions are 

V k (y) = Q, Vk= 1 , 2 , n ( 3 ) 

Where V k (y) = a k y(a ) + a^y\a) + a[ 2) y"{a ) + + a[ n ~ l) y u '~ l) (a) + J3 k y(b ) + p\y\b) 

+ J 8l 2) y\b)+ +j8 < k ”- ]> y ( "- ]) (b) (4) 

Where the linear forms Vj, V 2 . V n in 


y(a),y (a) . y {n n (a),y(b), y (b) . y ,n l \b) are linearly independent. 

Suppose the homogeneous boundary value problem given by (1) to (4) has only a trivial solution 
y(.r) = 0. The Green function of the boundary value problem (1) to (4) is the function G(x,t ) 
constructed for any point tin a<t <b and has the following properties 

1- In each of the interval [a, t) & (t, Z?]the function G(x,t )is considered as a function of x 

and is a solution of (1) i.e L(G) — 0. (5) 

2- G(x,t ) is continuous and has continuous partial derivatives with respect to x up to (n - 2) 
orders on a<x<b. 

3- The (n — 1) th derivative of G(x,t) with respect to x at x — t has a discontinuity of 1st kind 

f 1 'l 

and the jump being equal to-i.e 

l Po(OJ 

f d'-'G ') f 1 " 

dx"~ l dx"- 1 p 0 (O, 

V / x=t +0 V / x=t-0 V ruv 7 / 

4- G(x,t) satisfies the boundary conditions V k {G) = 0, k = 1,2.... n. 

Note- If the boundary value problem ( 1) to (4) has only a trivial solution y(x) = 0, then the 
operator L has a unique Green function G(x, t). 
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Example- Find the Green function of the following boundary value problems 


l- y=o, 


y(0) = y(/) = 0 


2 - y"+ju 2 y = 0, 


y(0) = y(l) = 0 




3- xy + y - — fy = 0, 

U/ 


y(0) is finite & y(l) = 0 


Solution 1- Let y =0 y(0) = y(l ) = 0 

First, we show that (1) has a zero or trivial solution. 

For this integrating (1) twice with respect to x, then the solution of (1) is 
y-ax+b 

Using the boundary conditions y(0) = 0 => b — 0 and y(l ) = 0=>0 = a/=>a=0 
Hence y(x) = Ois a zero or trivial solution of (1). So (1) has a Green function 


G(x,t) = ■ 


ja 1 x + a 2 0 <x<t 
[b l x + b 2 t<x<l 


Now G(x,t ) satisfies the following properties 

(1) G(x,t) is continuous at x — t so b l t + b 2 =a { t + a 2 => (b { -a x )t = {a 2 -b 2 ) 

1 


(2) G(x,t ) has a discontinuity of magnitude 
r dG^ fdG A 


V dx y x=t+0 


V & y x=f _ 0 


P 0 (0 

= — 1 => — a : = -1 


i.e 


(3) G(x,t ) satisfies the boundary conditions G(0,0 = 0 => a,.0 + a 2 =0=>a, =0 
G(/, t) = 0 => bj + b 2 =0 
From (3) & (4) a 2 -b 2 = —t 

Solving (3), (4), (5), (6) & (7) we have b 2 = t ,b, = ■-}, «, = 1 •-j 


Therefore G(x,r) = 


x(l—t) 


l 

t(l-x ) 


/ 


0 < x < r 
t < X<1 


, _ „ , . , sinr/r sin//(t-l) „ , sinwrcos// 

Answer 2- Solution a, = 0, Z?, = —— , a 2 =- — - &b 2 = -—-— 




p sin ju 


p sin // 



G(x,t) = • 


sin//(/-l) sin fix 
//sin // 

sin///sin//(/-l) 
//sin // 


0 < x < t 


t < x < 1 


Answer 3 Let xy + y - y = 0 


y(0) = finite y(l) = 0 


Solution- (1) can be written as x 2 y + xy -y= 0 


y(0) = finite = k(say) y(l) = 0 (2) 


i i 1 2 

Putting x = e" & x— = — = D, x 2 —j = D(D - 1) in (2), we have 
dx du dx~ 


(D 2 -l)y = 0 


The solution of (3) is 


y=ax+— 
x 


G(x, t ) = 


xY \-t 


2A t J 
Y l - x” 
2j[ x 


0 <x<t 


t < x < 1 


Sturm’s Comparison theorem: Let us consider two system of 2 nd order differential equations 
(pfit)u ) +qfit)u = 0 (1) & (p 2 (t)u) +q 2 (t)u = 0 (2) and assume that 

1- p l (t),p 2 (t), qfit) & g 2 (/)are continuous on the [a,b\. 

2- ufit) & u 2 (t ) are the non trivial solutions of the given system of equations. 

3- equation (2) is a sturm majorant of (1) on \a,b\. 

4- the inequality ^ > Pi( a ) u i( a ) holds, an q the LHS & RHS become infinite 

ufia) u 2 (a) 

when Mj (a) = 0 & u 2 (a) = 0. 

5- ufit) has exactly /7>lzerosat t = t l ,t 2 . t n where t 1 <t 2 <.</ H of ( a > b\. 

Then the solution M 2 (/)has at least n-zero on (a, t n \. 


Proof: Let us define a pair of continuous functions say (ftfit) & (f) 2 {t)by 
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<p { (t) = tan 


</> 2 (t) = tan 


»i (0 

Pi(0«i(0 y 
u 2 (t) 


p 2 (t)u 2 (t ) 

on the [a,Z?]for 0 < (Z \ <7t, V i = 1, 2 (only principle values). 


( 1 ) 

( 2 ) 


Since p. (?) > p 2 (?) =>-<-=>- ; — <- T -so for a particular value at ? = a, 

PiO) p 2 (0 P 2 (t)u 2 (?) 

we have — — ^ — hence 0 < (a) < (j) 2 {a) < n 
p l (a)u l (a) p 2 {a)u 2 {a ) 

Also by prufer’ s transformation 


(?) = — 1 7 “cos 2 <j) x + ^(?)sin 2 ^ 
Pi(0 


0 2 (?) = —-—cos 2 (f) 2 + g 2 (?)sin 1 (j) 1 
p 2 (t) 


(3) 

(4) 


Putting/.(?, 0) =—-—cos 2 (Zf +<7 (?)sin 2 V? = l,2 using in (3) & (4), we have 

/M'> 

$(0 = /,(b <f>) & (/> 2 (?) = / 2 (?, 0 ) => = f t (?, 0 ) v 1 = 1,2 (5) 

The solution of (3) is (?) with condition (a, ^ l (n , ))and the solution of (4) is (j) 2 (t) with 
condition (a, <p 2 (a)). 

Since />,(?) > p 2 (?) =>—-—<—-— V?e [a, b] so by (5), we have 
Pi(t) p 2 (t) 

{ ( t ) \ = f\(t, (/>)} <{(f) { = f 2 {t, (/))} on integrating, we have (p x (t) < (p 2 (t) V?e[a, Z?]so putting 
? = t n we have (p x (t n )<<p 2 (t n ) => nit < (j) 2 (t n ) since (j) x (?) has exactly n zeros on (a, b]. This implies 
u 2 (?) has at least n zeros on (a, t n ] (where t n =b). 


Corollary: The zeros of two linearly independent solutions u x (?) & u 2 (?) of 

(pj(?)?< ) +q l (t)u = 0 interlace, i.e between any two consecutive zeros of one solution there lies a 
zero of the other solution. 



(1)then 


Proof: Letter) & w 2 (t)be any two solutions of (p, (t)u ) +q l (t)u = 0 
both m, (t) & u 2 (t) satisfy ( 1 ) so 

{Pi (t)u'i ( 0 ) + q t ( t)u x (t) = 0 ( 2 ) 

{p x it)u 2 (t)) + q l (t)u 2 (t) = 0 (3) 

Let t 1 & t 2 be any two consecutive zeros of u 1 (t) so u l (t 1 ) = 0 = u 1 (t 2 ). 

Subtracting on multiplying (2) by u 2 and (3) by u t and then integrating between the limits to t 2 
, we have 

\p l (t){u' l (t)u 2 (t) - u 1 (t)u 1 (t)]\ r2 = 0 
Or 

L Pi (t 2 ){w, (t 2 )u 2 (t 2 ) - u 2 (t 2 )u x (t 2 )}J-1 p l (t i ){u[ (t, )u 2 o, ) - u 2 (t x )Mj (r, )}J = 0 
putting Mj (t { ) = 0 = M| (t 2 ), so it become 

Pl (t 2 )i/l (? 2 (^2Pi (?i )Wj (?i )t/2 (^l) • 

Since m^) & «j(f 2 ) are of opposite signs so u 2 (t ] ) & w 2 (t 2 ) will be of opposite signs. Hence 
m 2 (0 has at least one zero between t ] & t 2 . 



